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EXECUTIVE  StJMMARY 


Between  1982  and  1986,  the  Biogeochemical  Studies  program  of 
the  Acidic  Precipitation  in  Ontario  Study  included  in-depth 
investigations  of  the  interactions  between  incident  precipitation 
and  forest  trees  in  two  catchments  in  the  Muskoka-Haliburton  area 
of  Ontario.  These  investigations  were  conducted  under  a  contract 
with  the  Faculty  of  Forestry,  University  of  Toronto. 

The  catchments  had  contrasting  cover  characteristics,  and 
were  under  hydrologie  and  chemical  budget  monitoring  by  other 
studies.  Incident  precipitation,  throughfall  and  stemflow  were 
measured  and  sampled  by  distributing  collectors  within  the  basins 
using  stratified  and  fixed  plot  sampling  strategies.  Concentra- 
tions of  the  major  ions  in  the  samples  were  determined. 

The  forest  cover  characteristics  of  distinct  stands  on  the 
two  basins  are  described  in  terms  of  species  composition,  basal 
area,  stem  frequency  and  areal  extent. 

The  hydrologie  and  chemical  data  were  scrutinized  and  edited 
to  ensure  that  all  possible  observations  were  accounted  for,  con- 
formed to  chemical  laws  and  intuitive  principles,  and  agreed  with 
concomitant  observations.  This  report  describes  these  validation 
and  editorial  procedures. 

The  appendices  to  this  report  include  data  tables,  in  micro- 
computer file  format,  containing  hydrologie  and  chemical  charac- 
teristics of  incident  precipitation,  throughfall  and  stemflow 
collected  during  the  course  of  this  study. 
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1  INTRODUCTION 

In  1982,  the  Air  Resources  Branch  of  the  Ministry  of  the  ■ 
Environment  contracted  with  the  Faculty  of  Forestry,  University 
of  Toronto  to  undertake  investigations  on  two  forested  watersheds 
in  the  Muskoka-Haliburton  area  of  Ontario.  These  investigations 
were  conducted  under  the  auspices  of  the  Acidic  Precipitation  in 
Ontario  Study  (APIOS)  and  comprised  part  of  the  ' Biogeochemical 
Studies'  program  of  APIOS. 

A  major  component  of  these  investigations  was  the  monitoring 
of  hydrological  pathways  and  the  flux  of  materials  in  and  out  of 
the  watershed  as  well  as  between  watershed  compartments.  Such 
monitoring  included  the  collection  and  chemical  characterization 
of  precipitation  water  at  three  points  along  the  path  from  the 
atmosphere  to  the  forest  floor.  The  first  of  these  was  unaltered 
incident  precipitation.  The  second  point  was  throughfall,  defined 
as  precipitation  water  that  had  been  at  least  partially  inter- 
cepted by  the  forest  canopy  before  reaching  the  forest  floor.  The 
third  was  stemflow,  consisting  of  water  that  had  been  intercepted 
by  the  canopy  and  had  flowed  down  the  tree  bole. 

This  report  describes  the  methodologies  employed  in  monitor- 
ing incident  precipitation,  throughfall  and  stemflow  and  the 
rationale  behind  these  methods.  It  describes  the  techniques  used 
to  validate  data  and  generate  data  estimates  in  cases  were  the 
original  data  were  rejected  as  invalid  or  were  not  collected. 
Appendices  to  this  report  consist  of  MS-DOS  diskettes  containing 
data  tables. 

2  STUDY  SITES 

2.1  SELECTION  CRITERIA 

The  watershed  basins  selected  for  this  study  were  located  in 
the  Muskoka-Haliburton  area  of  south  central  Ontario  and  are 
designated  as  Plastic  Lake  No.  1  and  Harp  Lake  No.  4.  These 
basins  met  the  selection  criteria  which  were: 

•  The  basins  had  defined  catchment  areas  and  were  headwater 
basins  receiving  chemical  inputs  solely  from  atmospheric 
sources . 

•  The  basins  were  in  geographic  proximity  and  subjected  to  sim- 
ilar and  concurrent  météorologie  activity  and  atmospheric 
loading  of  long-range  transport  air  pollutants. 

•  The  basins  had  contrasting  characteristics  of  forest  cover, 
soil  type  and  depth  and  demonstrated  differences  in  the  water 
quality  of  streams  draining  the  basins. 
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•  The  streams  draining  the  basins  were  gauged  and  the  water 
quality  monitored. 

•  The  basins  had  reliable  all-season  access  and  reasonable 
assurance  of  tenure. 

2.2  PLASTIC  LAKE  NO.  1 

Plastic  Lake  is  located  at  latitude  45'  11'  N  and  longitude 
78'  49'  W  in  Sherborne  Township,  Haliburton  County.  The  lake 
surface  is  at  an  elevation  of  378  metres  a.s.l.   It  is  a  small 
Precambrian  Shield  headwater  lake,  32.14  hectares  in  size  and 
receives  drainage  from  a  92.50  hectare  catchment.  The  Plastic 
Lake  No.  1  basin  occupies  a  westerly  sector  of  the  lake  catch- 
ment and  is  23.34  hectares  in  size.  This  basin  is  drained  by  a 
perennial  stream  which  is  gauged  by  a  V-notch  weir  at  381 
metres.   The  maximum  elevation  in  the  basin  is  42  0  metres,  gen- 
erating a  relief  of  39  metres. 

The  basin  is  defined  primarily  by  natural  topography.  The 
exception  is  an  easterly  section  of  boundary  which  follows  a 
gravel  road  bed.  This  road  provides  vehicle  access  to  the  basin 
during  the  summer  and  fall  seasons.  Winter  access  is  by  snowmo- 
bile and  spring  access  is  by  all  terrain  cycle.  All  land  com- 
prising the  basin  is  managed  by  the  Leslie  M.  Frost  Natural 
Resources  Centre  of  the  Ontario  Ministry  of  Natural  Resources. 
There  is  no  development  on  any  part  of  the  lake  catchment.. 

The  bedrock  of  the  basin  is  dominated  by  gneiss  of  igneous 
origin  with  outcrops  of  metasedimentary  features.  Rock  outcrops 
are  frequent,  otherwise  the  surficial  geology  consists  of  thin 
basal  tills  on  the  uplands  and  organic  deposits  in  the  lowlands. 
The  majority  of  the  soils  belong  to  the  podzolic  order  and  vary 
from  zero  to  approximately  two  meters  in  depth. 

2.3  HARP  LAKE  NO.  4 

Harp  Lake  is  located  at  latitude  45'  23'  N  and  longitude  79' 
8 '  W  in  Chaffey  Township,  Muskoka  District  Municipality.  The 
lake  surface  is  at  an  elevation  of  327  metres  a.s.l.  It  is  also 
a  headwater  lake,  but  is  substantially  larger  at  71 . 38 ' hectares 
and  drains  a  470.66  hectare  catchment.  The  Harp  Lake  No.  4  basin 
is  119.09  hectares  in  size  and  occupies  a  westerly  sector  of  the 
lake  catchment.  The  basin  is  drained  by  a  perennial  stream  orig- 
inating in  a  pond  but  receiving  input  from  tributaries  along  its 
course  to  the  lake.  The  stream  is  gauged  by  a  rectangular  weir 
at  an  elevation  of  329  metres.  The  maximum  elevation  on  the 
basin  is  420  metres,  generating  a  relief  of  91  metres. 
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This  basin  is  also  primarily  defined  by  natural  topography 
except  for  a  relatively  short  section  of  gravel  road  bed  on  the 
easterly  extreme  of  the  basin.  This  road  can  be  used  by  conven- 
tional motor  vehicles  year  round.  All  land  comprising  the  basin 
is  privately  owned.  Harp  Lake  is  developed  with  cottages  and 
year-round  homes.  Some  sections  of  the  Harp  4  basin  have  been 
logged  and  are  un- forested. 

The  bedrock  of  this  basin  is  primarily  amphibolite  and 
schist.  The  surficial  geology  consists  of  basal  tills  on  the 
uplands  and  sand  deposits,  frequently  covered  by  peat,  in  the 
lowlands.  Bedrock  is  only  exposed  at  bedrock  ridges.  Podzolic 
soils  dominate  the  uplands,  with  organic  soils  occurring  in  the 
lowland  peat  areas.  Till  depth  is  highly  variable  and  can  exceed 
50  metres  in  valley  bottoms. 

3  FOREST  COVER 

3.1  SURVEY  TECHNIQUES 

Details  on  the  techniques  used  to  suirvey  the  standing  crop 
characteristics  of  the  Plastic  Lake  No.  1  and  Harp  Lake  No.  4 
basins  are  given  in  a  report  by  Lozano  and  Parton  (1987)^.  All 
forest  cover  data  presented  in  this  section  are  taken  or  adapted 
from  this  report. 

Briefly,  the  survey  consisted  of  delineating  a  fixed  grid 
system  on  the  watershed  and  using  a  wedge  prism  with  a  basal 
area  factor  of  2.0  m^/ha  to  determine  the  basal  area  and  stem 
density  of  different  species  at  each  grid  intersection.  This  was 
followed  by  caliper  measurements  of  each  stem  within  the  prism 
plot  for  stems  greater  than  10  cm  Diameter  at  Breast  Height 
(DBH)  at  Plastic  1  and  greater  than  12  cm  DBH  at  Harp  4  to 
provide  a  diameter  class  distribution  for  each  species. 

The  basins  were  also  divided  into  forest  cover  stands  based 
on  distinguishable  differences  in  species  composition  as  deter- 
mined by  air  photo  interpretation  and  ground  observations.  The 
Plastic  1  basin  was  divided  into  five  forest  cover  types  plus 
rock  outcrop  and  alder  swamp  areas.  The  Harp  4  basin  was  divided 
into  13  forest  cover  types  plus  un-forested,  open  water  and  road 
areas . 


1  Unpublished  Report:  Forest  Cover  Characteristics  of  the  Harp  #4 
and  Plastic  #1  Subcatchments  of  the  Southern  Ontario  Biogeochemi- 
cal  Study,  Faculty  of  Forestry,  University  of  Toronto 
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The  locations  of  each  of  the  stand  types  and  other  geo- . 
graphic  features  comprising  Plastic  Lake  No.  1  and  Harp  Lake  No. 
4  basins  are  given  in  Figures  1  and  2  respectively.  Also 
represented  in  these  figures  are  the  locations  on  the  basins  of 
the  incident  precipitation,  throughfall  and  stemflow  sampling 
stations  to  be  discussed  later. 

3.2  CHARACTERISTICS 

The  characteristics  of  each  forest  stand  type  are  summarized 
in  Tables  1  through  4 .  Two  approaches  are  used  in  these  summa- 
ries. In  Tables  1  and  2  the  basal  areas  occupied  by  the  various 
species  are  expressed  as  percentages  of  the  total  basal  area  of 
the  standing  crop.  In  Tables  3  and  4  the  frequencies  of  occur- 
rence of  individual  stems  of  each  species  are  summarized.  In 
both  approaches,  the  summaries  are  limited  by  exclusion  of 
individual  trees  smaller  than  the  minimum  DBHs  discussed  above. 

Tables  1  through  4  also  contain  summaries  of  the  character- 
istics of  the  total  forested  area  of  each  basin.  These  were 
derived  by  weighting  the  individual  stand  characteristics  by  the 
areas  occupied  by  those  stands.  Un-forested  land,  open  water, 
swamps,  bedrock  outcrops  and  roads  were  not  included  in  the 
weighting.  Areas  occupied  by  these  features  are  reported  sepa- 
rately. 

The  order  of  dominance  of  the  species  differs  depending  on 
whether  basal  area  percentage  or  stem  frequency  is  used  for 
ranking.  This  reflects  the  stem  diameter  distribution  for  the 
different  species.  Small  stem  diameter  species,  or  a  predomi- 
nance of  younger  trees  will  contribute  to  higher  frequency  val- 
ues, while  predominance  of  larger  species  or  older  trees  will 
contribute  to  higher  percent  basal  area  values. 

Basal  area  can  be  considered  a  better  indicator  of  the  domi- 
nance by  a  given  species  since  it  considers  both  the  number  and 
size  of  the  trees  and  intuitively  the  dominance  of  the  canopy. 

The  dominant  species  at  Plastic  1  are  White  Pine  and  Eastern 
Hemlock,  both  conifers.  Soft  Maple,  Red  Oak  and  White  Birch  are 
the  major  .deciduous  species  present.  White  Cedar  is  also  a  major 
component  of  the  basin  but  it  is  restricted  to  a  wet,  swampy, 
inaccessible  topographic  depression  and  to  the  gullies  leading 
to  this  swamp  area. 
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One  single  stand  type  constitutes  54.5  %  of  the  total  basin 
area.  Stand  A  is  dominated  by  the  two  major  conifer  species  with 
only  Soft  Maple  and  Red  Oak  as  significant  contributors  to  the 
stand.  White  Pine  and  Eastern  Hemlock  are  also  major  contrib- 
utors to  the  remaining  forested  stands  with  the  exception  of  the 
cedar  swamp  designated  as  Stand  C. 

Harp  4  has  been  divided  into  13  forest  stand  types.  The  size 
of  this  basin,  about  five  times  the  area  of  Plastic  1,  is  the 
major  reason  for  this  apparent  finer  stratification.  Most  of  the 
stand  types,  specifically  the  larger  ones,  are  dominated  by 
deciduous  species.  Hard  Maple  and  Soft  Maple  constitute  51.3  % 
of  the  basin  basal  area.  Largetooth  Aspen,  American  Beech  and 
Yellow  Birch  are  also  significant  components.  Eastern  Hemlock  is 
the  major  component  of  Stand  H  in  the  southwestern  part  of  the 
basin. 

The  differences  between  the  two  basins  can  be  summarized  as 
follows.  Plastic  Lake  No.  1  is  a  small  catchment  basin  charac- 
terized by  very  thin  podzolic  soils  and  dominated  by  White  Pine 
and  Eastern  Hemlock,  both  conifer  species.  Secondary  species  are 
various  hardwoods.  Harp  Lake  No.  4  is  a  much  larger  basin  with 
deeper  podzolic  soil  and  mixed  hardwoods,  which  are  dominated  by 
Hard  Maple  and  Soft  Maple.  Other  deciduous  species  contribute 
most  of  the  balance  of  the  standing  crop.  The  only  conifer  of 
significance  is  Eastern  Hemlock,  but  it  is  confined  to  a  unique 
forest  stand. 
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TABLE    1:    Plastic   Lake  No.    1    Forest   Cover   -   Basal   Area 


Stand 

Cover   Characteristics* 

Area    (ha) 

%   of   Total 

A 

PW47   Hei9   Or,   MSg   BW5   Oio 

12.72 

54.50 

B 

MS28   PW21   Cei2   Bwii   Org   Heg   0^0 

3.32 

14.22 

C 

Ce,i   Sbi3   Swg   Oio 

2.18 

9.34 

D 

PW32   MS23   Oris   Bwii   O16 

1.90 

8.14 

E 

He35   Ori4   Pwi4   Bw^^  Ms^   POg   O3 

1.80 

7.71 

F 

Exposed  Bedrock 

0.85 

3.64 

G 

Alder   Swamp 

0.57 

2.44 

Watershed  Basin  Totals 

23.34 

100.00 

PL    1** 

PW37  Hei7  MSi3   Or^o   Ce7   BW7   O9 

21.92 

93.92 

Subscript  number   following  the   species   code  represents   the  percentage 
of   the  basal   area   of   the   stand  or   forested  portion  of  watershed  con- 
tributed  to  by   that   species.     'Other'    species    include   those   contribut- 
ing less   than  five  percent   of   the  basal   area. 


Forested   stand   area  weighted  average    cover   characteristics. 


Species  Codes 

Bw 
Ce 
He 
Ms 

White   Birch 
White   Cedar 
Hemlock 
Soft  Maple 

Or 
Po 
Pw 

Red  Oak 
Poplar 
White  Pine 

Sb 
Sw 
0 

Black   Spruce 
White   Spruce 
Other   Species 
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TABLE  2:  Harp  Lake  No.  4   Forest  Cover  -  Basal  Area 


Stand 

Cover   Characteristics' 

Area    (ha) 

%   of   Total 

A 

Mh32   MSi8   Bei4   He^   O22 

21.78 

18.29 

B 

Mh48    P09   MSg    Bwdg   By-7    Fbs   0^5 

15.82 

13.28 

C 

Mh34   MSi6   PO14   Bwio   Beg   O^g 

14.33 

12.03 

D 

Mhjg   By2i   MSii    POg   BWj   O^e 

11.40 

9.57 

E 

Mh3-7   MS22   PW9    P07   Chbs   Be5   0^5 

9.31 

7.82 

F 

Mh^i  MS18   Byio   POio   Beg   0^6 

7.81 

6.56 

G 

Non-Forested 

6.79 

5.70 

H 

He26   MSi3    Sbi3   Fbn   Bwdio   Ce^o   O^, 

6.63 

5.57 

I 

Mhso   MSi7    Byg    Po,    0^^ 

5.14 

4.32 

J 

Mh82   Be5   Oi3 

4.55 

3.82 

K 

MS32   Mh29    Beii    Chbg   Bys   O^g 

3.47 

2.91 

L 

Mh53    Bwdii   By^i   AWg   Oig 

2.87 

2.41 

M 

Fb23   Ce2o    Swi4   Abio   By^o   Bwd,    BW5   0^ 

2.77 

2.33 

N 

PO44  Mhjo   Bei5   MSii   Oio 

1.38 

1.16 

0 

Surface  Water   &   Road 

5.04 

4.23 

Watershed  Basin  Totals 

119.09 

100.00 

HL    4*' 

Mhj,  Msi4   P07   Be,  By,  Heg  O23 

107.26 

90.07 

Subscript  number  following  the  species  code  represents  the  percentage 
of  the  basal  area  of  the  stand  or  forested  portion  of  watershed  con- 
tributed to  by  that  species.  'Other'  species  include  those  contribut- 
ing less  than  five  percent  of  the  basal  area. 


Forested  stand  area  weighted  average  cover  characteristics. 


Species  Codes 

Ab 

Black  Ash 

Ce 

Cedar 

Po 

Poplar 

Aw 

White  Ash 

Chb 

Black  Cherry 

Pw 

White   Pine 

Be 

Beech 

Fb 

Balsam  Fir 

Sb 

Black   Spruce 

Bw 

White  Birch 

He 

Hemlock 

Sw 

White   Spruce 

Bwd 

Basswood 

Mh 

Hard  Maple 

0 

Other   species 

By 

Yellow  Birch 

Ms 

Soft  Maple 
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TABLE   3:    Plastic   Lake  No.    1   Forest   Cover   -    Stem  Frequency 


Stand 

Cover   Characteristics* 

Area    (ha) 

%   of   Total 

A 

Pw29   He27   MSi7   Ce,   Org   BW5   Oio 

12.72 

54.50 

B 

MS37   Cei9   Pwii   Bwio   He,   Fbg   O9 

3.32 

14.22 

C 

Ce^s   Sbi2   SWg   Og 

2.18 

9.34 

D 

MS40   Pwi9   Bwio  Org   Heg    Iw^   0^ 

1.90 

8.14 

E 

He48  MSi5   Bwi4   Or9   PWg   0, 

1.80 

7.71 

F 

Exposed  Bedrock 

0.85 

3.64 

G 

Alder   Swamp 

0.57 

2.44 

Watershed  Basin  Totals 

23.34 

100.00 

PL    1" 

He22  PW21  MS20   Cei4   Bw,   Org   0^ 

21.92 

93.92 

Subscript  number   following   the  species   code  represents   the  stem  fre- 
quency of   that   species   in  the   stand  or   forested  portion  of   the 
watershed.     'Other'    species    include   those   contributing   less    than   five 
percent   of    the   stem   frequency. 


**       Forested  stand  area  weighted  average   cover   characteristics. 


Species  Codes 

Bw 
Ce 
He 
Ms 

White  Birch 
White  Cedar 
Hemlock 
Soft  Maple 

Or 
Po 
Pw 

Red  Oak 
Poplar 
White  Pine 

Sb 
Sw 
0 

Black   Spruce 
White   Spruce 
Other   Species 
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TABLE   4:    Harp   Lake  No.    4    Forest   Cover   -    Stem  Frequency 


Stand 

Cover  Characteristics* 

Area    (ha) 

%   of   Total 

A 

Mh33   MS22   Ben   Heg   Bwdg   O21 

21.78 

18.29 

B 

Mh52   MSg    POg    Fbg   Bwde   Bys   0^5 

15.82 

13.28 

C 

Mh35   MSi8    BWio    PO9    Beg    IWg    Oi4 

14.33 

12.03 

D 

Mh43    Byi9   MSg   Abs   O24 

11.40 

9.57 

E 

Mh35   MSi5    PWi2    POii   Chb,    BWg   O13 

9.31 

7.82 

F 

Mh^g  MSi4   Byg   Be,   P07   Abg   0^0 

7.81 

6.56 

G 

Non-Forested 

6.79 

5.70 

H 

He22   Fbi7   Bwdis  MSi4   Ce^   Sbg   0^3 

6.63 

5.57 

I 

Mh59   MSi3   By,    POg   Oi5 

5.14 

4.32 

J 

Mhgi    IW7   Be5   0, 

4.55 

3.82 

K 

Mh37   MS25   Beg   Chbg   Byg   O15 

3.47 

2.91 

L 

Mh53    Bwdii   Byio  AWg   O^j 

2.87 

2.41 

M 

Fb24    Ce25   Abi3    SWii   Bwdg   By,   BWg   0, 

2.77 

2.33 

N 

PO30   Mh24    Be23   MSi6   0, 

1.38 

1.16 

0 

Surface  Water   &   Road 

5.04 

4.23 

Watershed  Basin  Totals 

119.09 

100.00 

HL    4** 

Mhjg   MSi4    Bee    P05   Bys    O31 

107.26 

90.07 

Subscript  number   following   the   species   code  represents   the   stem  fre- 
quency of   that   species   in   the   stand  or   forested  portion  of   the 
watershed.    'Other'    species    include   those  contributing   less   than   five 
percent   of   the   stem  frequency. 


Forested  stand  area  weighted  average   cover  characteristics. 


Species   Codes 

Ab 

Black  Ash 

Ce 

Cedar 

Po 

Poplar 

Aw 

White  Ash 

Chb 

Black  Cherry- 

Pw 

White    Pine 

Be 

Beech 

Fb 

Balsam  Fir 

Sb 

Black   Spruce 

Bw 

White  Birch 

He 

Hemlock 

Sw 

White   Spruce 

Bwd 

Basswood 

Mh 

Hard  Maple 

0 

Other   species 

By 

Yellow  Birch 

Ms 

Soft  Maple 
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FIGURE  1:  Plastic  Lake  No.  1  Forest  Stand 
and  Collector  Locations 
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FIGURE  2:  Harp  Lake  No.  4  Forest  Stand 
and  Collector  Locations 
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4  COLLECTOR  DESIGN  AND  SAMPLING  STRATEGIES 

4.1  INCIDENT  PRECIPITATION 

Incident  precipitation  was  collected  by  two  types  of  collec- 
tors, depending  on  whether  the  anticipated  precipitation  was 
rain  or  snow.  Both  were  open  continuously  and  are  considered 
'bulk'  collectors. 

The  rain  collector  consisted  of  a  polyethylene  funnel,  28.8 
cm  in  diameter  finished  to  knife  edge  and  with  a  catch  area  of 
650  cm^.  The  funnel  stem  opening  was  covered  with  a  100  micron 
nylon  mesh.  The  stem  was  inserted  into  a  hole  in  the  lid  of  a 
five  litre  polyethylene  cylinder  and  sealed  with  laboratory 
film.  The  inside  of  the  cylinder  was  lined  with  a  bag  consisting 
of  a  polyethylene  interior,  laminated  to  a  nylon  exterior  to 
improve  puncture  resistance.  Two  identical  collectors  were  used 
in  the  collection  of  incident  rain. 

The  snow  collector  consisted  of  a  80  litre,  polyethylene 
refuse  container  lined  with  a  larger  polyethylene /nylon  bag.  A 
polyethylene  ring  was  sized  to  fit  the  inside  perimeter  of  the 
top  of  the  container.  This  ring  kept  the  bag  in  place  and  pro- 
vided a  catch  area  of  1600  cm^ . 

Incident  precipitation  was  collected  above  the  canopy  by 
erecting  a  tower  of  construction  scaffolding  and  placing  the 
collectors  on  the  top  platform.  At  Plastic  1,  the  tower  was 
located  at  a  bedrock  outcrop,  while  at  Harp  4  the  tower  was 
located  on  a  small  hilltop.  The  collectors  were  at  an  elevation 
of  15  meters  above  the  ground.  Figures  1  and  2  indicate  the 
locations  of  these  towers  in  the  respective  basins. 

Construction  of  the  towers  permitted  the  placement  of  inci- 
dent collectors  in  proximity  to  the  respective  throughfall  col- 
lectors without  relying  on  or  creating  a  clear  area.  The 
contamination  of  incident  precipitation  samples  by  wind  blown 
tree  litter  (leaves,  seeds,  pollen)  was  also  reduced. 

In  addition  to  the  two  incident  rain  funnels,  a  standard 
rain  gauge  and  a  battery-operated,  tipping  rain  gauge  were  also 
installed  on  the  towers.  The  standard  gauge  would  be  used  to 
resolve  any  funnel  collected  depth  inconsistencies,  while  the 
tipping  bucket  would  provide  a  continuous  record  of  rain  activ- 
ity; start  and  end  times,  intensity,  duration.  During  the  snow 
season,  incident  precipitation  was  collected  by  a  single  80 
litre  snow  collector.  The  gauges  were  not  appropriate  for  winter 
use. 
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At  the  Plastic  1  bedrock  outcrop,  a  standard  rain  gauge  or 
one  of  the  snow  collectors  was  also  installed  at  ground  level 
near  the  tower.  The  surrounding  forest  provided  a  much-needed 
wind  screen,  especially  for  the  snow  collector.  The  precipita- 
tion depths  measured  by  this  collector  would  eventually  displace 
the  measurements  recorded  by  the  snow  collector  on  the  tower. 
The  latter,  it  was  determined,  underestimated  the  amount  of 
incident  snow. 

A  similar  bedrock  outcrop  or  other  type  of  forest  clearing 
was  not  available  at  the  Harp  4  basin. 

Throughout  the  rain  seasons,  incident  precipitation  was 
sampled  by  'event'.  An  event  was  defined  as  a  period  of  precipi- 
tation activity  isolated  by  24  hours  of  dry  weather.  During  the 
snow  season,  collections  were  made  weekly,  given  that  there  was 
an  accumulation  of  precipitation  during  the  preceding  week  and 
there  was  no  precipitation  occurring  on  the  scheduled  collection 
day.  Otherwise  the  collection  was  made  the  following  week  or 
when  activity  ceased. 

The  collection  consisted  of  removing  and  sealing  the  liner 
bags.  On-site  weighing  of  the  sample  was  performed  with  a  spring 
scale  to  provide  a  more  reliable  volume  estimate  in  case  of  sam- 
ple loss  during  transport.  New  bags  were  installed  and  either 
the  funnel  or  snow  collector  ring  replaced.  At  10  to  14  day 
intervals,  lab-washed  rain  funnels  were  installed  during  the 
collection  visits  to  prevent  the  accumulation  of  organic  materi- 
als, tars  and  resins,  from  the  trees  on  the  funnel  surfaces. 
Washing  was  performed  with  laboratory  detergent  followed  by 
several  rinses  with  deionized  water. 

4.2  THROU6HFALL 

The  throughfall  collectors  were  identical  to  the  incident 
precipitation  collectors  in  design.  Also,  sample  collection, 
collector  maintenance  and  volume  recording  procedures  were  iden- 
tical and  were  performed  concurrently  with  those  for  the  inci- 
dent precipitation. 

During  the  course  of  this  study,  two  distinct  throughfall 
sample  collector  placement  strategies  were  employed.  The  first 
strategy  involved  positioning  the  collectors  beneath  'target' 
canopies  that  represented  the  range  of  species  and  diameter 
classes  in  the  basin.  This  stratified  approach  had  two  phases, 
an  intense  phase  and  a  reduced  phase  where  the  number  of 
throughfall  collection  locations  was  lower.  The  second  strategy 
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utilized  a  fixed  plot  in  the  different,  major  stand  types.  In 
addition,  within  the  reduced  stratified  sampling  period,  sam- 
pling beneath  replicated  target  canopies  was  employed. 

The  aim  of  each  strategy  was  to  characterize  the  role  of  the 
whole  of  the  forested  area  of  each  basin  in  altering  both  the 
quantity  and  quality  of  precipitation  reaching  the  forest  floor. 
The  exception  to  this  objective  occurred  with  the  exclusion  of 
the  cedar  swamp  area.  Stand  H,  at  Plastic  1.  This  stand  was  too 
inaccessible  for  routine  monitoring. 

The  first  stage  in  designing  the  stratified  approach  con- 
sisted of  determining  the  number  of  throughfall  collector  loca- 
tions which  could  be  attended  to  with  finite  logistic  resources. 
Twenty  five  stations  at  each  basin  were  the  maximum  which  could 
be  established.  The  second  phase  was  to  determine  the  canopy 
types  which  could  be  addressed  with  these  collectors.  At  the 
Plastic  1  basin,  the  four  dominant  species  were  selected  as  tar- 
get canopy  types.  At  Harp  4,  the  five  most  dominant  species  were 
selected.  Dominance  was  defined  as  higher  basal  area  for  reasons 
discussed  above.  These  four  and  five  species  collectively  com- 
prised 77  and  72  percent  of  the  basal  area  of  the  forest  trees 
in  the  Plastic  1  and  Harp  4  basins  respectively.  Finally  the 
collector  locations  were  stratified  in  proportions  that  repre- 
sented the  relative  dominance  of  each  of  the  selected  species  on 
the  watershed  and  encompassed  the  range  of  stem  diameters 
exhibited  by  that  species. 

Two  or  three  (for  larger  trees)  of  the  f\innel-type  rain  col- 
lectors were  positioned  within  the  drip  line  of  each  target 
canopy.  The  volume  of  each  individual  collector  was  determined, 
however  the  samples  would  be  bulked  prior  to  chemical  analysis. 
A  single  snow  throughfall  collector  was  similarly  positioned. 
The  position  of  each  independent  collector  was  marked  so  that  it 
could  be  returned  to  the  same  position  at  the  beginning  of  the 
given  season  if  required. 

Tables  5  and  6  detail  the  distribution  of  throughfall  col- 
lectors during  the  different  phases  of  the  stratified  strategy. 
The  alphabetic  codes  refers  to  the  target  species  and  the 
sequential  numbers  identifying  each  station  within  the  species 
group.  Within  each  group,  the  numbers  represent  increasingly 
larger  diameters.  Stations  identified  with  the  characters  'b' 
and  'c'  represent  the  replicates  of  certain  target  canopies 
which  were  established  to  measure  the  degree  of  variability  in 
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throughfall  depth  and  chemical  composition  that  may  occur  under 
target  canopies  of  like  species  and  similar  stem  diameter  and 
crown  shape  and  size. 

This  sampling  strategy  should  not  be  construed  as  attempting 
to  isolate  the  role  of  a  given  tree  species.  The  influence  of 
the  overstory  at  locations  targeting  smaller  diameter  trees  or 
the  influence  of  the  understory  near  the  larger  trees  was  pres- 
ent and  no  attempt  was  made  to  reduce  it.  The  strategy  attempted 
to  ensure  that  collectors  were  placed  in  a  variety  of 
micro-cover  types.  This  would  permit  a  simple  extrapolation  of 
the  observations  made  at  the  selected  locations  to  the  whole 
basin,  with  the  exception  of  certain  stands  previously  dis- 
cussed. 

The  stratified  sampling  was  suspended  and  replaced  by  a 
fixed  plot  strategy  in  April,  1985.  Six  plot  locations  were 
established  in  each  basin.  At  Plastic  1,  two  plots  were  estab- 
lished in  each  of  stands  A,  D  and  E.  At  Harp  4,  one  plot  was 
established  in  each  of  stands  A,  B,  D  and  F.  Two  plots  were 
established  in  stand  C.  Figures  1  and  2  show  the  locations  of 
these  plots  in  the  respective  basins. 

The  Plastic  1  stands  selected  to  contain  the  plots,  collec- 
tively resembled  the  overall  cover  characteristics  of  the  basin. 
On  an  areal  basis,  stands  A,  D  and  E  accounted  for  70  %  of  the 
basin  and  75  %  of  the  forested  area.  The  two  stands  excluded 
from  this  coverage,  while  representing  the  second  and  third 
largest  stands,  were  unique  and  presented  logistic  problems. 
Both  of  these  stands  occupied  low  lying  areas,  the  gully  and  the 
adjacent  slopes,  and  the  cedar  swamp. 

The  Harp  4  stands  A,  B,  C,  D  and  F  likewise  represented  the 
overall  cover  of  the  basin  and  included  the  largest  stands. 
These  stands  represented  60  %  of  the  basin  area  and  66  %  of  the 
forested  area.  Stand  E  was  excluded  for  access  reasons.  Stand  C 
received  two  plots  because  it  in  itself  best  represented  the 
basin  forest  cover.  Harp  4  also  contained  a  number  of  minor  yet 
distinct  stands  which  could  not  be  addressed  with  finite 
resources.  Most  of  these  were  hardwood  dominated  as  were  the 
primary  stands.  Their  distiction  lay  in  the  relative  dominance 
of  different  species.  The  exceptions  were  two  stands,  H  and  M, 
with  conifer  dominance. 
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The  location  of  each  plot  within  a  stand  was  selected  to 
represent  a  'typical'  cover  condition  in  that  stand  with  respect 
to  species  composition,  tree  size  and  stocking  density.  The 
plots  were  15  x  15  meters  and  were  staked  at  3 -meter  intervals. 
Five  rain  or  three  snow  collectors  were  randomly  re-distributed 
among  the  stake  points  during  each  clean  rain  funnel  exchange  or 
during  each  snow  period  collection.  The  volume  of  each  individ- 
ual collector  was  determined,  however  the  samples  were  bulked 
prior  to  any  chemical  analysis. 

•  Tables  5  and  6  contain  dates  with  which  the  first  and  last 
collection  in  each  strategy  or  phase  was  identified.  A  Reference 
Number  identifies  each  collection  over  the  term  of  the  study. 
This  designation  will  be  discussed  later.  Finally,  these  tables 
report  the  number  of  collections  made,  the  number  of  throughfall 
stations  active  and  the  total  number  of  samples  collected. 

4.3  STEMFLOW 

The  design  of  the  stemflow  collector  consisted  of  a  length 
of  laboratory  tubing,  brand  name  'Tygon  Formulation  2506',  of  1" 
O.D.  and  3/4"  I.D.  This  tubing  was  split  longitudinally,  wrapped 
in  a  downward  spiral  around  the  tree  bole,  fastened  with  stain- 
less steel  staples  and  sealed  to  the  bark  with  acrylic  caulking. 
On  species  with  rough  bark,  the  bark  was  shaved  to  improve  the 
seal  while  taking  care  not  to  damage  the  cambium  layer.  The 
lower,  unsplit  end  of  the  tubing  was  inserted  into  a  hole  in  the 
lid  of  a  plastic  refuse  container  which  was  lined  with  a  dispos- 
able polyethylene  refuse  bag.  The  capacity  of  the  container  was 
125  litres. 

The  collection  of  samples  involved  weighing  the  liner  bag 
with  a  spring  scale,  or  if  the  sample  was  too  large,  using  a 
calibrated  plastic  rod  to  measure  the  volume  of  the  sample.  The 
sample  was  stirred  to  ensure  complete  mixing  and  an  aliquot  was 
taken  into  a  475  mL  polystyrene  jar.  Excess  sample  was  poured 
onto  the  ground  at  the  base  of  the  respective  bole. 

The  strategies  for  sampling  of  stemflow  paralleled  those  of 
throughfall.  There  was,  however  was  no  sampling  of  stemflow  dur- 
ing the  fixed  plot  throughfall  sampling  period.  Stemflow  collec- 
tors were  used  on  the  bole  of  each  of  stratified  period  target 
canopies,  including  the  replicated  canopies,  at  the  same  time 
and  for  the  same  duration. 
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To  reduce  the  logistic  demands  on  the  field  and  laboratory- 
staff,  a  restriction  on  the  definition  of  a  stemflow  sample  was 
placed.  It  was  assumed  that  a  sample  existed  if  more  than  one 
litre  was  collected.  The  rational  behind  this  restriction  was 
that  the  bole  beneath  the  collar  would  have  absorbed  some  of  the 
stemflow  water  if  it  had  not  been  intercepted.  All  stemflow  vol- 
umes were  reduced  by  this  amount  to  maintain  consistency.  This 
restriction  was  admittedly  arbitrary,  but  in  light  of  the  large 
volumes  of  such  samples,  it  was  unlikely  to  present  a  noticeable 
bias  to  the  data.  It  did,  however,  ensure  that  each  stemflow 
sample  would  be  of  sufficient  size  for  complete  chemical  analy- 
sis . 

5  SAMPLE  PROCESSING  AND  CHEMICAL  ANALYSIS 

Upon  collection  and  weighing,  all  samples  were  transported 
to  the  field  laboratory  at  the  Dorset  Research  Center  and  pro- 
cessed for  chemical  analysis. 

The  weight  of  the  contents  of  each  bag  was  determined  on  an 
electronic  pan  balance  to  improve  accuracy  of  this  measure  prior 
to  bulking.  In  cases  of  large  quantities  of  sample,  exceeding  the 
capacity  of  the  lab  balance,  the  weight  determined  in  the  field 
was  used.  Aliquots  of  incident  rain  and  throughfall  were  trans- 
ferred to  disposable  475  mL  polystyrene  jars. 

These  jars,  along  with  those  containing  the  stemflow  were 
delivered  to  the  analytical  laboratory  at  the  research  center 
where  the  pH  and  Gran  acidity  were  determined.  These  measurements 
were  made  within  the  same  day  to  within  five  days  of  sample  col- 
lection in  the  field  depending  on  logistic  constraints.  The 
remaining  sample  was  shipped  to  the  Laboratory  Services  Branch 
(LSB)  laboratory  in  Toronto  for  determination  of  major  ions. 

Details  pertaining  to  the  analytical  methods  can  be  found  in 
various  analytical  methodology  dociiments  developed  by  LSB.  In 
brief,  pH  was  determined  with  a  pH  meter  and  combination  glass- 
calomel  electrode  and  Gran  acidity  by  titration  with  a  base  to  a 
fixed  end-point.  The  anions  sulphate,  nitrate  and  chloride  were 
determined  by  ion  chromatography,  while  the  base  cations  calcium, 
magnesium,  potassium  and  sodiiom  were  determined  by  flame  atomic 
absorbtion  spectrophotometry.  Ammonixim  ions  were  determined  by  a 
colourimetric  technique.  Phosphate  ions  were  also  determined  but 
their  concentrations  were  too  frequently  below  detection  limits 
to  be  useful. 
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TABLE  5:  Sampling  Periods  and  Strategies  at  Plastic  Lake  No.  1 


SAMPLE  DATE  (first) 

830105 

840215 

840418 

840903 

850419 

SAMPLE  DATE  (last) 

840208 

840411 

840830 

850416 

860415 

REF.  NO.  (first) 

36 

104 

112 

146 

187 

REF.  NO.  (last) 

103 

111 

144 

185 

255 

Incident 

PW  1 

PW  2 

PW  3 

PW  3b 

PW  3c 

PW  4 

PW  5 

PW  6 

PW  7 

PW  7b 

PW  7c 

PW  8 

PW  9 

PW  10 

PW  11 

PW  lib 

PW  lie 

PW  12 

HE  1 

HE  2 

HE  3 

HE  3b 

HE  3c 

HE  4 

HE  5 

HE  6 

HE  6b 

HE  6c 

MS  1 

MS  2 

MS  2b 

MS  3 

MS  4 

MS  4b 

OR  1 

OR  2 

OR  3 

Plot  1 

* 

Plot  2 

* 

Plot  3 

•. 

Plot  4 

* 

Plot  5 

* 

Plot  6 

* 

Collections 

64 

8 

28 

35 

61 

Incident  Samples 

64 

8 

28 

35 

61 

Throughfall  Stations 

25 

8 

20 

8 

6 

Throughfall  Samples 

1600 

64 

560 

280 

366 

Stemflow  Stations 

25 

8 

20 

8 

0 

Stemflow  Samples 

556 

31 

246 

124 

0 

Stemflow  not  collected 
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TABLE  6:  Sampling  Periods  and  Strategies  at  Harp  Lake  No.  4 


SAMPLE  DATE  (first) 

830105 

840215 

840418 

840831 

850428 

SAMPLE  DATE  (last) 

840208 

840411 

840830 

850423 

860424 

REF.  NO.  (first) 

36 

104 

112 

145 

188 

REF.  NO.  (last) 

103 

111 

144 

186 

256 

Incident 

MH  1 

MH  lb 

MH  Ic 

MH  2 

MH  3 

MH  4 

MH  5 

MH  5b 

MH  5c 

MH  6 

MH  7 

MH  8 

MH  9 

MH  10 

MH  11 

MH  lib 

MH  lie 

MH  12 

MH  13 

* 

MS  1 

MS  2 

MS  2b 

MS  3 

MS  4 

MS  4b 

MS  5 

MS  6 

PO  1 

PO  2 

BY  1 

BY  2 

BE  1 

BE  2 

Plot  1 

• 

Plot  2 

• 

Plot  3 

* 

Plot  4 

* 

Plot  5 

* 

Plot  6 

* 

Collections 

66 

8 

26 

40 

62 

Incident  Samples 

66 

8 

26 

40 

62 

Throughfall  Stations 

25 

8 

16 

8 

6 

Throughfall  Samples 

1650 

64 

416 

320 

372 

Stemflow  Stations 

24 

8 

16 

8 

0 

Stemflow  Samples 

814 

39 

275 

178 

0 

■*  =  stemflow  not  collected 
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6  THE  DATA  BASE 

6.1  SIZE 

Tables  5  and  6,  which  describe  the  collection  strategies 
employed  in  this  study,  also  report  the  nimber  of  collections 
and  active  stations  used  to  monitor  incident  precipitation, 
throughfall  and  stemflow,  at  each  catchment  basin.  Also  reported 
are  the  number  of  samples  resulting  from  these  collections.  This 
information  is  reported  for  each  collection  strategy. 

Note  that  the  number  of  incident  and  throughfall  samples  is 
a  product  of  the  number  of  stations  and  collections.  This 
assumes  that  throughfall  occurred  at  each  station  during  every 
sampling  period  and  that  there  were  no  accidental  losses  of  sam- 
ple. There  were  in  fact  rare  cases  when  there  were  no  through- 
fall  samples  as  a  result  of  a  very  light  precipitation  event.  A 
'sample'  of  depth  equal  to  zero  represents  such  cases  in  the 
data  base.  Where  a  sample  was  accidently  lost,  a  record  was 
created  to  generate  such  a  sample.  The  procedure  for  dealing 
with  accidental  losses  is  discussed  in  the  next  section. 

The  number  of  stemflow  samples  represents  the  actual  number 
of  cases  were  stemflow  of  more  than  one  litre  occurred.  The  num- 
ber of  locations  (boles)  with  stemflow  could  vary  from  zero  to 
the  maximum  of  stations  active  during  the  collection  period, 
depending  on  the  characteristics  of  the  precipitation.  The  data 
base  size  statistics  include  cases  where  estimates  of  stemflow 
volume  and  chemistry  were  created  to  compensate  for  the  acciden- 
tal sample  loss.  The  procedures  used  to  recreate  a  lost  sample 
are  also  discussed  in  the  next  section.  Note  that  stemflow  was 
not  sampled  during  the  fixed  plot  strategy.  Also,  Station  MH  13 
at  Harp  4  did  contain  a  stemflow  collar,  however  the  tree  bole 
contained  a  major  crack  and  the  collar  could  not  be  attached 
properly.  This  resulted  in  consistently  unreliable  stemflow  vol- 
ume data.  All  stemflow  data  collected  at  this  station  were 
removed  from  the  working  database. 

6.2  ORGANIZATION 

Two  key  fields  are  required  to  identify  any  sample  collected 
during  the  course  of  this  study.  These  are  the  Sample  Date  and 
Location  Code.  During  the  periods  employing  the  rain  funnels, 
the  Sample  Date  represents  the  calendar  day  during  which  the 
last  part  of  the  'event'  occurred.  It  does  not  represent  the  day 
the  samples  were  removed  from  the  field.  During  the  period 
employing  the  snow  collectors,  the  Sample  Date  does  represent 
the  day  the  samples  were  removed  from  the  field.  Recall  that 
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monitoring  of  precipitation  characteristics  by  the  tipping 
bucket  gauge  and  recorder  was  not  possible  during  the  winter  and 
it  was  not  possible  to  record  the  time  when  precipitation 
stopped.  All  samples,  incident,  throughfall  and  stemflow,  col- 
lected concurrently,  have  the  same  Sample  Date. 

The  second  key  field  is  the  Location  Code.  This  field  can' 
best  be  explained  by  the  following  examples. 

'PR-INC   represents  incident  precipitation  collected  by  the 
rain  funnel  at  the  Plastic  Lake  basin. 

'HS-MS4b'  represents  throughfall  collected  by  the  snow  collec- 
tor under  the  target  canopy  Soft  Maple  Number  4b  at 
the  Harp  Lake  Basin.  The  character  'b'  identifies 
this  as  one  of  the  replicated  canopies. 

'PS-PT3'   represents  throughfall  collected  by  snow  collectors 
in  Plot  3  at  the  Plastic  Lake  basin. 

'HSF-BE2'  represents  stemflow  collected  from  target  canopy 
American  Beech  number  2  at  the  Harp  Lake  basin. 

To  facilitate  the  recognition  of  each  sampling  period,  all 
sample  collections  were  assigned  a  sequential  collection  Refer- 
ence Number  extending  over  the  entire  study  period.  The  same 
n\imber  was  allocated  to  each  sample  collected  at  the  respective 
basin  at  any  time.  The  same  number  occurring  for  collections  at 
both  basins  identifies  a  concurrent  period  of  precipitation 
activity.  Collections  with  the  same  Reference  Number  may  have 
■different  Sample  Dates.  The  most  frequent  cause  of  this  apparent 
discrepancy  is  the  cessation  of  precipitation  on  either  side  of 
midnight  at  the  two  basins. 

A  number  applying  to  only  one  of  the  basins,  signifies  that 
precipitation  activity  did  not  occur  at  the  other  basin  or  that 
the  samples  do  not  represent  concurrent  activity  or  collection 
period.  Some  numbers  do  not  apply  to  either  basin.  This  signi- 
fies that  precipitation  was  collected  only  at  the  Dorset 
Research  Centre  where  an  incident  collector  was  co-located  with 
collectors  of  other  studies  for  intercomparison  purposes.  This 
referencing  was  used  mainly  in  the  validation  procedures  and 
data  management . 

The  Reference  Numbers  begin  at  number  36.  Monitoring  of  some 
nature  began  at  Plastic  1  in  early  August,  1982,  and  in  late 
August,  1982  at  Harp  4.  Since  this  was  a  new  study,  it  was 
necessary  to  develop  procedures  and  correct  problems  as  the  they 
became  apparent.  Consequently,  data  collected  during  the  early 

SDB-016-3512-93  21 


phases  of  collection  were  of  questionable  quality.  Once,  collec- 
tion with  snow  buckets  began,  a  new  set  of  problems  had  to  be 
resolved.  It  was  decided  to  exclude  all  collections  with  a 
Reference  Number  less  than  3  6  from  the  data  base,  and  therefore 
from  any  data  analysis  or  interpretation. 

6.3  PRELIMINARY  EDITING 

During  the  course  of  the  study,  data  was  continuously  writ- 
ten to  mainframe  computer  files.  These  data  consisted  of  sample 
dates,  numeric  sampling  location  codes,  sample  numbers,  sample 
volumes  and  a  series  of  chemical  concentrations  for  each  sample. 
These  files  were  continuously  edited  to  ensure  that  the  correct 
field  and  field  laboratory  information  had  been  transcribed. 
This  information  included  the  sample  date,  location  code,  sample 
number,  sample  volume  and  pH.  Field  and  field  lab  notes  were 
compared  to  the  contents  of  the  mainframe  files  to  ensure  con- 
sistency. The  balance  of  the  chemical  data,  which  was  generated 
by  the  main  laboratory,  was  also  reviewed  to  identify  apparent 
errors.  If  a  transcription  or  laboratory  error  was  confirmed, 
appropriate  corrections  were  made  to  the  mainframe  files. 

In  addition,  where  no  information  was  available  for  a  sample 
due  to  accidental  loss,  an  artificial  record  for  this  sample  was 
created  to  ensure  that  it  was  not  overlooked.  Estimates  for  the 
missing  data  would  be  generated  later. 

Manual  screening  of  mainframe  data  file  reports  identified 
chemical  data  that  were  clearly  inconsistent  with  other  data 
from  similar,  concurrent  samples.  These  gross  outliers  could 
most  easily  be  spotted  in  a  listing  of  all  throughfall  samples 
collected  at  one  time.  Such  suspect  data  could  be  a  result  of 
sample  contamination  or  lab  instrumentation  and  reporting 
errors.  These  were  deleted  from  the  mainframe  files.  Also 
deleted  were  chemical  data  reported  as  being  greater  than  some 
value,  since  such  information  had  limited  value,  and  would 
create  a  bias  if  left  in  the  data  base  as  a  real  concentration. 
Estimates  for  these  deleted  data  would  also  be  created  later. 

7  DATA  BASE  EDITING  AND  VALIDATION 

7.1  RATIONALE 

As  a  result  of  such  editing,  the  mainframe  files  contained 
records  for  each  sample  collected,  or  that  should  have  been  col- 
lected, and  these  records  were  correctly  identified.  However, 
many  of  these  records  were  incomplete,  with  some  or  all 
observations  missing,  and  otherwise  not  completely  validated. 
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All  available  mainframe  data  were  transferred  to  microcomputer 
applications  for  further  editing,  validation  and  creation  of 
estimates  for  the  voids. 

The  occurrence  of  erroneous  or  missing  data  could  be 
ascribed  to  several  causes.  These  are  discussed  below  in 
decreasing  order  of  frequency. 

The  most  frequent  cause  of  gaps  in  the  data  was  due  to 
insufficient  sample  volume  for  a  complete  chemical  analysis. 
Since  sampling  was  conducted  on  an  event  basis,  it  was  not 
appropriate  to  wait  for  the  subsequent  precipitation  activity  to 
add  sufficient  volume  to  the  samples. 

Certain  samples,  especially  those  of  small  volume,  could 
contain  concentrations  of  analytes  which  exceeded  the  analytical 
range  of  the  instrumentation.  If  there  was  insufficient  sample 
for  dilution  and  re-analysis  the  datum  was  reported  as  being 
greater  than  the  instrumentation  limit.  These  data  had  been 
deleted  from  the  data  base  as  mentioned  above. 

The  third  most  common  cause  of  missing  information  was  the 
accidental  losses  of  samples.  This  could  occur  before  or  during 
collection  resulting  in  no  information  about  the  sample.  A  loss 
could  occur  during  transport  to  the  field  lab,  also  resulting  in 
no  information  unless  a  gravimetric  determination  of  sample  vol- 
ume had  been  made  at  collection.  Samples  could  also  be  lost  at 
any  stage  of  the  chemical  analysis  resulting  in  only  partial 
chemical  information. 

Recording  and  transcription  errors,  either  by  the  field  of 
laboratory  personnel,  contamination  of  or  improper  handling  or 
preservation  of  samples,  or  laboratory  instrument  problems  would 
result  in  erroneous  information  in  the  data  base.  Such  errors 
were  not  always  apparent  during  the  preliminary  screening  and 
edit . 

A  significant  source  of  error  was  related  to  the  poor  catch 
efficiency  of  the  snow  collectors  which  were  situated  on  the 
towers.  This  highly  exposed  environment  was  subject  to  wind  tur- 
bulence above  the  collector  resulting  in  the  capture  of  less 
than  the  actual  incident  depth  of  snow. 

Finally,  there  was  the  possibility  of  contamination  of  a 
sample  during  collection  by  field  crews,  or  in  the  transfer  to 
various  containers.  Foreign  substances  could  also  be  present  in 
the  sample.  This  would  be  especially  evident  in  incident  precip- 
itation samples. 
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7.2  PROCEDURES 

7.2.1  INCIDENT  PRECIPITATION  DEPTH 

Since  the  data  on  depth  of  incident  precipitation  occurring 
at  the  study  sites  were  critical  to  a  valid  interpretation  of 
the  throughfall  process,  it  was  vital  that  all  erroneous  or 
missing  information  for  this  parameter  be  replaced  by  the  best 
estimate  possible. 

During  the  rain  seasons,  precipitation  was  collected  by 
paired  funnels  on  the  tower.  If  the  depths  of  the  funnels, 
which  were  averaged,  were  not  consistent  or  did  not  agree  with 
the  standard  gauge,  the  record  was  edited  to  reflect  the  higher 
funnel  depth.  This  procedure  would  have  to  improve  the  agree- 
ment between  funnel  and  gauge  depths.  If  both  funnel  depths 
were  not  consistent  with  the  standard  gauge  the  later  was  used. 
The  remark  '  dbl  '  ,  representing  a  doubling  of  the  higher  fiinnel 
depth  qualifies  the  depth  datum.  The  remark  'std'  represents 
cases  were  the  standard  gauge  datum  was  used.  Further  evidence 
of  an  erroneous  incident  rain  depth  was  the  occurrence  of  sev- 
eral throughfall  depths  which  were  higher  than  the  respective 
incident  depth,  an  intuitively  unacceptable  condition. 

Such  changes  also  usually  resulted  in  better  agreement  with 
precipitation  depth  observations  made  by  an  independent  study 
operating  in  the  basins,  as  well  as  eliminating  the  throughfall 
versus  incident  depth  inconsistency. 

This  independent  study,  which  shall  be  referred  to  as  the 
WRB  study,  maintained  incident  collectors  at  the  two  basins  as 
part  of  an  investigation  of  hydrological  and  chemical  budgets. 
The  sampling  periods  used  by  the  WRB  study  did  not  coincide 
with  those  of  this  study,  and  were  invariably  longer  periods. 
Consequently,  two  or  more  of  either  study  sampling  periods  were 
combined  until  a  coincident  sampling  date  was  reached.  This  was 
not  a  difficult  task  since  both  studies  collected  samples  after 
the  cessation  of  precipitation  and  common  sampling  dates  were 
frequent. 

During  the  snow  season,  there  was  a  recognized  problem  of 
poor  catch  efficiency  by  the  un-shielded  collectors  on  the  tow- 
ers. At  the  Plastic  basin,  a  separate  ground  level  snow  collec- 
tor was  located  on  the  un- forested  rock  outcrop  and  adjacent  to 
the  tower.  The  surrounding  forest  provided  an  effective  wind 
break  for  the  collection  of  snow.  Most  of  the  incident  snow 
depths  for  this  basin  were  based  on  this  collector  and  have  the 
remark  'grd'  appended.  This  substitution  invariably  resulted  in 
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a  good  agreement  with  the  independent  WRB  study  depth  records, 
while  at  the  same  time  eliminated  the  incidence  of  higher  than 
incident  throughfall  depths. 

Such  a  ground  level  collector  was  not  available  at  the  Harp 
basin  while  poor  catch  efficiency  was  also  a  problem.  The  inci- 
dent snow  depths  as  recorded  on  the  Hairp  tower  were  consis- 
tently lower  than  the  concurrent  throughfall  depths  and  lower 
than  the  WRB  study  data. 

To  compensate  for  these  erroneous  incident  snow  depth  data, 
histograms  of  the  incident  and  throughfall  depths  were  gener- 
ated for  each  collection  and  intuitive  estimates  of  valid 
depths  made.  While  this  procedure  was  being  used,  consideration 
was  given  to  the  type  of  snow  fall  that  had  occurred.  This 
included  the  total  amount,  consistency,  intensity  and  the  prob- 
ability that  under  some  conditions  the  snow  would  be  captured 
by  the  canopy  only  to  be  released  during  the  subsequent 
collection  period.  Invariably  when  the  winter  precipitation  was 
in  the  form  of  rain  or  wet  snow,  the  catch  efficiency  would  be 
good  and  the  throughfall  depths  would  be  consistently  less  than 
the  incident.  In  such  cases,  the  actual  collected  depth  was 
retained.  While  these  estimates  were  being  generated,  the  new 
depth  values  would  be  checked  for  agreement  with  the  WRB  study 
data.  The  remark  'tda'  was  appended  to  the  estimated  incident 
snow  depths  to  represent  a  throughfall  depth  adjusted  value. 

This  method  was  also  used  for  three  of  the  Plastic  basin 
incident  snow  depth  estimates  were  the  ground  level  collector 
was  disturbed  or  its  data  was  considered  unreliable. 

On  three  occasions,  the  WRB  study  depth  was  used  to  replace 
Harp  basin  incident  snow  depth  when  other  procedures  could  not 
be  applied.  In  these  cases  the  sampling  periods  of  the  two  stu- 
dies were  concurrent.  These  depth  data  have  been  identified 
with  the  remark  'wrb' . 

7.2.2  THROUGHFALL  DEPTH 

No  post   facto   changes  were  made  to  the  throughfall  rain  or 
snow  depth  data.  The  position  of  these  collectors  in  a  forest 
stand  provided  an  ideal  wind  break  so  that  the  collected  depths 
would  be  an  accurate  representation  of  below  canopy  precipita- 
tion penetration. 
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However,  over  the  course  of  the  study  there  were  cases 
where  a  given  throughfall  station  collector  was  disturbed  or 
depth  data  was  incomplete.  During  the  rain  periods,  this  con- 
sisted of  a  loss  of  information  from  one  or  more  of  the 
throughfall  funnels  positioned  under  a  target  canopy  or  from 
one  or  more  of  the  funnels  on  the  plot  grids.  The  depth  data 
was  then  solely  derived  from  the  remaining  collector (s) .  This 
reliance  on  a  single  collector  was  also  made  during  the  snow 
period  when  the  plot  strategy  was  used  since  there  were  three 
snow  collectors  used  at  each  plot.  The  remark  'E'  identifies 
all  throughfall  depth  data  that  was  obtained  in  a  less  than 
optimum  manner. 

In  a  few  cases,  there  was  no  throughfall  depth  data  for  a 
given  station  at  all.  This  occurred  primarily  where  only  a 
single  snow  throughfall  collector  was  used  beneath  a  target 
canopy  and  it  was  disturbed  or  the  depth  was  not  recorded.  In 
these  situations,  a  surrogate  collection  period  or  rain  event 
was  identified.  This  surrogate  would  have  similar  incident 
depth  and  would  be  as  close  as  possible  in  time  to  the  event 
with  the  missing  datum.  A  ratio  consisting  of  the  summed 
throughfall  depths  available  for  both  events  was  calculated  and 
this  ratio  was  applied  to  the  missing  datum  station  in  the  sur- 
rogate event  to  generate  a  value  for  the  missing  data  point. 
The  remark  'E'  was  also  used  to  identify  all  throughfall  depths 
so  derived. 

Using  the  above  described  procedures  resulted  in  a  complete 
and  more  accurate  record  of  incident  and  throughfall  depths  for 
all  stations  throughout  the  duration  of  the  study.  This  elimi- 
nated the  need  for  a  general  extrapolation  to  account  for  any 
missing  depth  data  which  would  not  have  considered  individual 
precipitation  characteristics  or  other  variables. 

7.2.3  STEMFLOW  VOLUME 

Correcting  erroneous  or  estimating  missing  stemflow  volumes 
presented  a  unique  problem.  Unlike  incident  or  throughfall 
depth  which  could  be  compared  to  independent  but  related  obser- 
vations, stemflow  volumes  could  only  be  obtained  from  actual 
observations.  If  a  sample  was  partially  or  completely  spilled, 
if  the  collar  was  blocked  by  ice  or  debris,  or  if  a  collector 
overflowed,  the  volume  had  to  be  estimated. 
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A  flagged  preliminary  estimate  was  added  to  the  mainframe 
files  as  a  place-holder  for  the  sample  record.  After  the  data 
were  transferred  to  a  microcomputer  environment,  all  such  esti- 
mates were  deleted.  Since  all  sampling  was  complete,  the  data 
base  contained  the  maximum  number  of  valid  observations. 

All  sampling  periods  that  contained  complete  and  valid 
stemflow  volumes  were  isolated  and  the  total  volume  of  stemflow 
for  each  bole  was  calculated.  To  produce  an  estimate  for  a 
missing  volume  within  a  sampling  period,  a  series  of  ratios 
between  the  individual,  valid  voliimes  of  that  sampling  period 
and  the  respective  total  voliomes  were  calculated.  The  average 
of  these  ratios  was  applied  to  total  stemflow  volume  for  that 
bole  to  generate  the  estimate  for  the  missing  stemflow  volume. 

In  this  way,  the  estimate  was  based  on  both  the  long  term 
stemflow  record  from  the  given  tree  bole  as  well  as  on  the 
flows  from  the  other  boles  for  the  given  sampling  period.  For 
this  approach  to  be  operational  required  grouping  of  sampling 
periods  having  the  same  sampling  strategy  or  phase.  The  proce- 
dure would  operate  as  long  as  there  was  at  least  one  valid 
volume  observation  for  each  sampling  period.  The  remark  'sfe' 
was  appended  to  all  stemflow  volumes  derived  in  this  way. 

7.2.4  INCIDENT  PRECIPITATION  CHEMISTRY 

Considerable  effort  was  expended  in  the  validation  or 
interpolation  of  incident  precipitation  chemical  data.  Two 
broad  approaches  were  used  in  these  procedures.  The  first 
involved  checks  of  the  internal  consistency  of  the  chemistry  of 
a  sample.  The  second  involved  comparisons  of  data  for  samples 
collected  by  other  collectors  and  by  other  studies,  at  nearby 
locations. 

One  internal  integrity  test  checked  the  ion  balance  of  a 
sample.  Since  the  analyses  performed  on  these  samples  included 
all  major  ions,  and  it  was  unlikely  that  incident  precipitation 
samples  contained  appreciable  quantities  of  dissociated  organic 
compounds,  it  was  possible  to  determine  the  discrepancy  between 
anion  and  cation  concentrations.  Consequently,  if  this  test 
indicated  that  the  anions  and  cations  were  in  balance,  then  it 
was  unlikely  that  any  one  concentration  was  in  error.  The  fol- 
lowing calculations  were  performed  to  determine  the  discrepancy 
between  the  reported  anion  and  cation  concentrations.  All 
concentrations  are  reported  in  milligrams  per  litre. 
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Anion  equivalence  -  Cations  e<iuivalence 

%   Discrepancy  =  : x  200 

Anion  equivalence  +  Cations  equivalence 

where: 

Anion  equivalence  = 

[SOJ       [NO3-N]        [CI]  4.90 


48.0288     14.0067      35.453     1000  x  10'"°** 

Cation  equivalence  = 

,  „„,3,      [NH4-N]       [Ca]         [K]         [Mg]        [Na] 


14.0067      20.04      39.0983     12.1525     22.9898 

The  second  internal  integrity  test  checked  the  measured  pH 
against  a  pH  value  that  was  predicted  from  the  concentrations 
of  other  ions  in  the  sample.  This  test  was  especially  valuable 
since  pH  was  potentially  the  least  stable  chemical  measurement 
made.  The  calculations  accommodated  dissolved  atmospheric  CO2 
in  the  sample . 

The  discrepancy  between  measured  pH  and  the  pH  predicted  by 
the  ionic  composition  of  the  incident  precipitation  sample  was 
determined  as  follows: 

%  Discrepancy  =  x  200 

P^aasurad  "•■  P^radictad 

The  predicted  pH  was  derived  from  the  calculated  hydrogen 
ion  concentration.  With  all  ion  concentrations  expressed  in 
milligrams  per  litre,  the  latter  was  determined  as  follows: 

(Net  Ions)  +  ({Net  lons)^  +  4(Kh  x  K,  x  PCO,))^'' 

[H+]  =   

2 

where  : 


Net  Ions  = 


[SOJ       [NO3-N]       [CI]       [NH,-N]       [Ca] 


48.0288     14.0067     35.453     14.0067      20.04 
[K]         [Mg]         [Na] 


39.0983     12.1525     22.9898 


X    1000 


and:  K^    =    0.034    X    10      ueq/L/atm       Henry's  Law  Constant 

Ki    =    0.45    ueq/L  First   CO2  Dissociation  Constant 

PCO2    =    320    X    10     atm  CO2  Partial   Pressure 

These   two   internal   integrity  tests  would  not   identify 
sample  contamination  by  an  ionic   compound  consisting  of  a  major 
anion  and  cation.    To   check   for   such  occurrences   and   to  produce 
estimates    for  missing  data,    data    for   concurrent  precipitation 
events    in   the   same   geographic   area  was   consulted. 
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As  mentioned  previously,  each  collection  period  for  the  two 
study  basins  was  assigned  a  reference  number.  An  identical  num- 
ber indicated  a  concurrent  sampling  period.  In  addition  to 
these  two  locations,  identical  collectors  were  operated  at  the 
Dorset  Research  Centre  (DRC)  by  this  study.  These  sampling 
periods  were  similarly  referenced,  a  relatively  easy  task  since 
all  three  stations  were  usually  sampled  on  the  same  day. 

A  fourth  set  of  data  was  available  from  an  automated  wet 
only  collector,  also  located  at  the  DRC  and  operated  by  the 
APIOS  Deposition  Network.  These  samples  were  collected  daily 
and  as  a  result  there  were  usually  more  than  one  sample  coin- 
ciding with  this  study's  sampling  period.  To  make  use  of  this 
fourth  data  set  required  the  identification  of  the  daily 
samples  which  coincided  with  the  sampling  periods  for  this 
study  and  calculating  the  depth  weighted  concentrations  for 
these  series  of  daily  samples.  The  resultant  groups  of  daily 
samples  were  assigned  the  same  series  of  reference  numbers. 

By  using  the  common  referencing  scheme,  it  was  possible  to 
compare  the  chemistry  of  any  one  incident  sample  to  three  oth- 
ers concurrently  sampled.  This  was  accomplished  by  merging  all 
four  data  sets  and  comparing  individual  data  from  each  set.  A 
concentration  reported  for  an  incident  precipitation  sample 
collected  at  one  of  the  basins  was  compared  to  concomitant  data 
from  the  other  sources.  If  precipitation  characteristics,  pri- 
marily depth,  were  similar,  the  suspect  datum  was  flagged. 

An  additional  source  of  information  for  validating  incident 
precipitation  chemistry  resided  in  the  chemistry  data  for 
throughfall  samples  collected  along  with  the  incident  sample  at 
the  basin.  Certain  assumptions  could  be  drawn  from  the  major 
trends  between  the  incident  and  throughfall  chemistry.  These 
were  that  sulphate,  chloride,  calcium,  magnesium,  potassium  and 
sodium  were  enriched  in  throughfall.  The  other  ions,  namely 
hydrogen,  nitrate  and  ammonium  could  be  enriched  or  reduced. 
Consequently,  if  the  reported  throughfall  concentrations  for 
the  ions  that  were  expected  to  be  enriched  were  consistently 
lower  than  in  the  incident  sample,  the  incident  sample  concen- 
tration was  considered  suspect. 

The  editing  and  validation  of  the  incident  precipitation 
chemistry  was  performed  by  using  all  the  tests  and  sources  of 
validating  information  simultaneously.  Electronic  spreadsheets 
containing  the  incident  data  from  each  basin  and  the  formulae 
for  internal  integrity  checks  were  prepared.  A  common  spread- 
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sheet  containing  merged  data  from  all  four  sources  was  likewise 
on  hand,  as  was  a  listing  of  incident  data  which  conflicted 
with  throughfall  data. 

The  data  files  contain  the  remark  "C"  identifying  incident 
precipitation  chemistry  data  which  were  changed  from  that 
reported  from  the  chemical  analysis.  The  remark  "E"  is  used  to 
identify  estimated  data  where  there  was  no  analytical  result 
reported. 

Data  installed  into  the  database  and  flagged  with  "C"  or 
"E"  had  to  be  compatible  with  all  validation  processes.  Thus, 
the  data  could  not  conflict  with  the  throughfall,  had  to  agree 
with  the  other  incident  precipitation  data  sources  and  hadi  to 
result  in  agreement  in  the  internal  integrity  checks. 

The  results  of  these  validation  and  editing  procedures  was 
to  have  a  chemical  concentration  value  for  each  major  ion  and 
each  sample  in  the  incident  precipitation  database  for  each 
basin. 

7.2.5  THROUGHFALL  CHEMISTRY 

All  throughfall  chemical  data  was  manually  screened  to 
identify  cases  where  values  were  outliers  from  the  range 
created  by  the  other  throughfall  locations.  This  was  the  second 
such  screening,  in  addition  to  the  preliminary  screening  con- 
ducted while  the  data  was  resident  on  the  mainframe  computer. 
The  remaining  outliers  were  deleted. 

The  remaining  throughfall  data  was  now  considered  valid 
although  there  where  a  large  number  of  voids  created  by  the 
deletions  and  where  no  data  were  generated.  Internal  integrity 
checks  were  not  performed  because  it  was  known  that  the  chemi- 
cal analyses  performed  were  not  inclusive  enough  to  cover  all 
possible  ions  in  such  samples.  Organic  anions,  for  example, 
would  have  a  significant  contribution  to  the  ion  balance  and 
predicted  pH  calculations.  These  were  not  analyzed. 

The  chemical  data  voids  were  filled  by  using  a  ratio  esti- 
mation technique.  This  technique  is  described  below. 

All  collection  periods  were  separated  on  the  basis  of 
throughfall  sampling  strategy  or  phase.  This  provided  data  sets 
consisting  of  the  same  series  of  collector  locations.  Further 
segregation  was  made  on  the  basis  of  whether  the  rain  or  snow 
collectors  were  used.  This  isolated  possible  seasonal 
influences . 
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Data  sets  within  each  new  grouping  which  did  not  contain 
any  voids  for  a  given  analyte  were  identified  and  simple,  un- 
weighted, grand  mean  concentrations  for  each  collector  location 
were  calculated. 

When  a  data  set  with  a  missing  value  was  encountered,  a 
series  of  ratios  between  available  concentrations  and  the 
respective  grand  mean  concentrations  were  generated.  The  aver- 
age of  these  ratios  was  then  calculated. 

This  average  ratio  was  then  applied  to  the  grand  mean  con- 
centration corresponding  to  the  collector  location  with  the 
data  void  to  produce  an  estimated  concentration. 

This  technique  had  the  advantage  of  taking  into  consider- 
ation the  chemistry  of  throughfall  collected  at  a  given  loca- 
tion over  the  course  of  a  collection  strategy  or  phase  as  well 
as  considering  the  chemistry  of  throughfall  collected  at  other 
locations  during  the  same  collection  period.  Throughfall  chem- 
istry concentrations  generated  by  this  technique  were  added  to 
the  data  base  and  flagged  with  the  remark  'RT'. 

To  use  this  technique  required  at  least  one  valid  chemical 
concentration  for  one  sampler  location  within  each  data  set. 
There  were  cases  where  all  concentration  data  for  a  certain 
chemical  were  missing.  This  occurred  when  very  small  volumes  of 
sample  were  collected  and  there  was  insufficient  sample  for 
analysis  of  all  chemicals.  This  required  a  different  approach 
to  produce  estimates. 

The  approach  adopted  consisted  of  conducting  linear  regres- 
sion analyses  between  all  ion  pairs  from  each  sampling  location 
within  the  sampling  strategy  or  phase.  The  ion  which  had  the 
best  correlation  with  the  missing  ion  was  used  in  the  regres- 
sion equation  to  calculate  the  missing  value.  Values  so  esti- 
mated were  installed  in  the  data  base  with  the  remark  'RG' . 

These  two  procedures  completely  filled  the  chemical  data 
base  for  major  ions  in  throughfall. 

7.2.6  STEMFLOW  CHEMISTRY 

Stemflow  chemistry  was  notable  in  the  wide  range  of  concen- 
trations which  could  occur.  A  preliminary  examination  of  the 
data  suggested  that  the  quality  of  stemflow  was  dependant  on 
total  flow  as  well  as  the  tree  species.  The  validation  and 
estimation  procedures  would  be  based  on  these  considerations. 
Procedures  such  as  the  ratio  or  regression  techniques  applied 
to  estimate  missing  throughfall  chemistry  would  be  inappropri- 
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ate  because,  unlike  throughfall  depths  which  fell  within  a  nar- 
row range  for  each  precipitation  period,  stemflow  volumes  could 
vary  by  orders  of  magnitude. 

The  procedure  used  to  fill  missing  stemflow  chemistry  data 
was  to  isolate  each  unique  bole  for  the  entire  study  period, 
without  regard  to  the  sampling  strategy.  A  correlation  matrix 
was  generated  for  each  bole  containing  the  correlation  coeffi- 
cients between  major  ions  plus  volume.  Whenever  a  stemflow 
record  with  a  missing  ion  concentration  was  encountered,  the 
matrix  was  examined  to  determine  which  other  ion  or  volume  had 
the  best  correlation  with  the  missing  ion.  If  the  stemflow 
record  contained  a  valid  concentration  for  this  other  ion,  the 
missing  ion  concentration  would  be  calculated  using  this  other 
ion  as  the  required  variable. 

Frequently,  the  other  ion  was  not  the  first  choice  based  on 
the  matrix  because  the  stemflow  record  had  several  missing  ion 
concentration  data.  At  times  there  were  no  ion  data,  only  a 
volume.  In  such  cases  all  missing  ion  concentrations  were  cal- 
culated using  volume  as  the  variable. 

If  the  correlation  between  the  ions  was  significant  (p  = 
0.05),  the  calculated  ion  concentration  was  added  to  the  data 
base.  This  procedure  was  repeated  to  fill  more  missing  data, 
but  on  a  more  complete  data  base.  The  remark  codes  'RGl'  to 
'RG6'  represent  the  first  through  sixth  passes  using  this  pro- 
cedure, at  which  each  datum  was  added.  Subsequent  passes  failed 
to  add  to  the  data  base.  This  left  a  relatively  small  number  of 
chemical  data  that  had  to  be  estimated  in  another  manner. 

These  data  were  generated  by  identifying  a  stemflow  sample 
from  the  same  tree,  with  a  similar  volume  and  occurring  during 
the  same  season.  The  appropriate  datum  was  copied  from  this 
sample  record  to  the  incomplete  sample  record.  The  remark  code 
'SUR'  identifies  such  copied  data. 

8  CALCULATION  OF  MEAN  CONCENTRATIONS 

The  data  files  appended  to  this  report  contain  ion  concen- 
tration data  for  the  various  samples  collected  during  the  course 
of  this  study.  To  reduce  these  data  to  manageable  proportions 
will  require  some  mathematical  manipulations.  Following  are  sug- 
gested approaches  for  data  reduction  based  on  the  sampling  design 
principles  of  this  study. 
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8.1  INCIDENT  PRECIPITATION 

Single  measurements  represent  depths  of  rain  'events'  or 
weekly  snowfall,  and  concentrations  of  major  ions.  Calculating 
mean  ion  concentrations  requires  weighting  the  individual  con- 
centrations by  the  respective  depths  for  the  period  of  interest. 

8.2  THR0U6HFALL 

The  throughfall  sampling  strategies  employed  in  this  study 
were  designed  to  simplify  the  assessment  of  the  influence  of  the 
watershed  forests  on  precipitation  quantity  and  quality.  The 
stratified  placement  of  throughfall  collectors  considered  basal 
area  dominance  of  tree  species.  Likewise,  the  placement  of  col- 
lectors in  fixed  plots  considered  the  stand  types  on  the 
watersheds . 

Calculation  of  mean  ion  concentrations  in  throughfall  for 
each  sampling  period  can  be  accomplished  by  depth  weighting  the 
individual  ion  concentrations  for  all  collector  locations.  This 
can  be  applied  to  the  intensive  and  reduced  stratified  sampling 
strategies.  It  should  be  recalled  that  during  the  reduced  stra- 
tified strategy,  there  was  a  period  where  certain  sampling  loca- 
tions were  'replicated' .  The  data  from  these  locations  should  be 
given  separate  consideration  and  not  included  in  depth  weighting 
calculations  unless  averages  for  these  replicates  are  first 
determined. 

The  throughfall  sampled  by  the  plot  strategy  requires  spe- 
cial consideration  if  the  total  watershed  forest  is  under  con- 
sideration. At  Plastic  Lake  No.  1,  six  plots  were  established, 
two  each  in  Stands  'A','D'  and  'E'.  At  Harp  Lake  No.  4,  two 
plots  were  established  in  Stand  'C  and  one  each  in  Stands  'A', 
'B',  'D'  and  'F'.  Since  the  stands  represent  varying  proportions 
of  the  watersheds,  throughfall  as  collected  at  the  different 
plots  should  be  weighted  when  determining  the  throughfall  char- 
acteristics of  the  portions  of  watershed  forest  under  investiga- 
tion. 

Table  7  presents  the  weighting  factors  to  be  applied  to  the 
plot  throughfall  data  to  improve  the  representation  of  the  whole 
watershed  forest.  Note,  however,  that  neither  the  whole  of  the 
watersheds  nor  their  forested  portions  are  completely  repre- 
sented under  any  sampling  scheme. 
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TABLE  7:  Weighting  Factors  for  Plot  Throughfall  Collections 


Harp  Lake  No.  4 

Plastic  Lake  No.  1 

Plot 
No. 

Stand 

Stand 
Area 
(ha) 

Weighting 
Factor 

Plot 
No. 

Stand 

Stand 
Area 
(ha) 

Weighting 
Factor 

PT  1 
PT  2 
PT  3 
PT  4 
PT  5 
PT  6 

B 
F 
C 
C 
D 
A 

15.82 
7.81 
14.33 
(14.33) 
11.40 
21.78 

0.2224 

0.1098 

0.1007 

0.1007  • 

0.1602 

0.3062 

PT  1 
PT  2 
PT  3 
PT  4 
PT  5 
PT  6 

D 
D 
E 
E 
A 
A 

(1.90) 
1.90 

(1.80) 
1.80 

12.72 

(12.72) 

0.0579 
0.0579 
0.0548 
0.0548 
0.3873 
0.3873 

Total 

71.14 

1.0000 

Total 

16.42 

1.0000 

8.3  STEMFLOW 

The  magnitude  of  flow  of  precipitation  water  down  the  boles 
of  trees  is  highly  dependant  on  the  bark  characteristics  (hence 
species)  and  size  of  the  tree  (DBH) .  Also,  different  species 
appear  to  have  greatly  differing  effects  on  the  chemistry  of  the 
water  that  flow  down  the  bole. 

Consequently,  the  calculation  of  mean  stemflow  concentra- 
tions as  well  as  determining  the  significance  of  stemflow  as  a 
hydrologie  pathway  requires  detailed  information  on  the  forest 
characteristics.  This  includes  species  composition,  size  class 
distribution  as  well  as  density  of  trees  in  various  species/size 
classes . 

Tables  8  and  9  summarize  the  distributions  of  tree  stems  in 
the  two  watershed  basins.  The  DBHs  of  the  sample  trees  are 
given,  along  with  the  range  of  DBHs  of  trees  of  that  species  in 
the  basin  meant  to  be  represented  by  the  sample  tree.  The  den- 
sity of  a  given  size  range  of  trees  in  stems  per  hectare  is 
required  information  to  calculate  fluxes  of  chemicals  via 
stemflow.  The  weighting  factor  is  a  ratio  of  the  density  of 
trees  in  the  size  class  to  the  density  of  all  trees  species  and 
size  classes  listed. 
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TABLE    8:    Tree   Stem  Distribution   -    Plastic   Lake  No. 


INTENSIVE  STRATIFIED 

STRATEGY 

Tree  Bole 

DBH  (cm) 

Stand  Range* 
DBH  (cm) 

S terns /Ha 

Weighting 
Factor 

PWl 

15.5 

10  -  16 

26.85 

0.04721 

PW2 

16.9 

18 

8.28 

0.01456 

PW3 

19.7 

20 

11.84 

0.02082 

PW4 

22.3 

22 

14.24 

0.02504 

PW5 

24.6 

24  -  26 

31.17 

0.05480 

PW6 

29.1 

28 

11.65 

0.02048 

PW7 

29.5 

30 

11.57 

0.02034 

PW8 

31.8 

32  -  34 

19.56 

0.03439 

PW9 

36.8 

36  -  38 

13.93 

0.02449 

PWIO 

39.5 

40 

5.93 

0.01043 

PWll 

43.9 

42  -  48 

16.53 

0.02906 

PWl  2 

52.7 

50  + 

3.09 

0.00543 

HEl 

13.6 

10  -  14 

71.75 

0.12614 

HE  2 

15.4 

16 

23.34 

0.04103 

HE  3 

17.4 

18 

16.71 

0.02938 

HE  4 

21.3 

20  -  22 

38.44 

0.06758 

HE5 

23.7 

24  -  26 

17.07 

0.03001 

HE  6 

29.3 

28  + 

18.61 

0.03272 

MSI 

13.8 

10  -  14 

63.83 

0.11222 

MS2 

18.1 

16  -  18 

56.21 

0.09882 

MSB 

21.4 

20  -  22 

26.53 

0.04664 

MS  4 

23.2 

24  + 

18.48 

0.03249 

ORl 

22.4 

10  -  26 

20.69 

0.03638 

0R2 

30.5 

28  -  34 

10.67 

0.01876 

0R3 

40.1 

36  + 

11.82 

0.02078 

REDUCED 

STRATIFIED  STRATEGY 

Tree  Bole 

DBH  (cm) 

Stanci  Range* 
DBH  (cm) 

St ems /Ha 

Weighting 
Factor 

PW3 

19.7 

10  -  22 

61.21 

0.10761 

PW7 

29.5 

24  -  34 

73.95 

0.13001 

PWll 

43.9 

36  + 

39.48 

0.06941 

HE3 

17.4 

10  -  22 

150.24 

0.26414 

HE  6 

29.3 

24  + 

35.68 

0.06273 

MS  2 

18.1 

10  -  18 

120.04 

0.21104 

MS4 

23.2 

20  + 

45.01 

0.07913 

0R2 

30.5 

10  + 

43.18 

0.07592 

DBH  range   of   trees   on   the  watershed  represented  by  designated  tree  bole 
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TABLE    9:    Tree   Stem  Distribution   -   Harp   Lake  No.    4 


INTENSIVE  STRATIFIED 

STRATEGY 

Tree  Bole 

DBH  (cm) 

Stand  Range* 
DBH  (cm) 

Stems/Ha 

Weighting 
Factor 

MHl 

13.3 

12  -  14 

71.85 

0.17955 

MH2 

15.3 

16 

27.91 

0.06975 

MH3 

18.5 

18 

26.23 

0.06555 

MH4 

19.8 

20 

23.02 

0.05753 

MH5 

22.0 

22 

16.57 

0.04141 

MH6 

24.3 

24 

15.33 

0.03831 

MH7 

26.3 

26 

10.41 

0.02601 

MH8 

28.0 

28 

8.41 

0.02102 

MH9 

30.5 

30 

6.50 

0.01624 

MHIO 

31.5 

32 

5.33 

0.01332 

MHll 

35.8 

34  -  38 

8.41 

0.02102 

MH12 

41.5 

40  + 

4.64 

0.01160 

MSI 

13.0 

12  -  16 

32.23 

0.08054 

MS2 

20.0 

18  -  20 

18.02 

0.04503 

MSB 

22.3 

22 

7.97 

0.01992 

MS4 

25.0 

24  -  26 

10.37 

0.02591 

MS  5 

28.0 

28  -  32 

7.97 

0.01992 

MS  6 

39.3 

34  + 

5.12 

0.01279 

POl 

24.3 

12  -  30 

26.62 

0.06652 

P02 

36.8 

32  + 

4.90 

0.01224 

BYl 

19.8 

12  -  26 

21.98 

0.05493 

BY2 

32.3 

28  + 

8.19 

0.02047 

BEI 

20.0 

12  -  26 

24.26 

0.06062 

BE2 

34.0 

28  + 

7.93 

0.01982 

REDUCED 

STRATIFIED  STRATEGY 

Tree  Bole 

DBH  (cm) 

Stand  Range* 
DBH  (cm) 

Stems /Ha 

Weighting 
Factor 

MHl 

13.3 

12  -  16 

99.76 

0.24929 

MH5 

22.0 

18  -  28 

99.97 

0.24982 

MHll 

35.8 

30  + 

24.88 

0.06217 

MS2 

20.0 

12  -  22 

58.22 

0.14549 

MS4 

25.0 

24  + 

23.46 

0.05863 

POl 

24.3 

12  + 

31.52 

0.07877 

BY2 

32.3 

12  + 

30.17 

0.07539 

BEI 

20.0 

12  + 

32.19 

0.08044 

DBH   range   of    trees   on   the  watershed  represented  by   designated   tree  bole 
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9  APPENDICES 

The  appendices  to  this  report  consist  of  eight  spreadsheet 
files  prepared  in  an  MS  DOS  microcomputer  environment  using  Lotus 
1-2-3,  ver.  2.x  software.  The  contents  of  each  file  are  described 
below.  They  are  configured  to  print  in  a  landscape  orientation  on 
a  Hewlett-Packard  Laserjet  series  printer.  The  files  are  con- 
tained on  three  3.5  ",  1.44  Mb  floppy  disks.  Sample  pages  from 
these  files  are  attached  to  this  report.  Copies  of  these  disks 
and  this  report  are  available  from: 

Ministry  of  Environment  and  Energy 

Dorset  Research  Centre 

Bellwood  Acres  Road 

Dorset,  Ontario 

POA  lEO 

Attn:  B.  LaZerte 

9.1  PIiASTIC    LAKE    NO.     1    INCIDENT   AND   THROUGHFALL 

Three  files  contain  incident  precipitation  and  throughfall 
data  collected  at  Plastic  Lake  No.  1  during  the  course  of  this 
study.  The  data  have  been  compiled  in  the  separate  files  along 
sampling  strategies.  File  names  and  periods  of  coverage  are  as 
follows: 

PL_TF_1.WK1       Intensive  Stratified  -  Reference  Nos.  36  -  103 
PL_TF_2.WK1      Reduced  Stratified  -  Reference  Nos.  104  -  185 
PL_TF_3.WK1       Plots  -  Reference  Nos.  187  -  255 

9.2  HARP  LAKE  NO.  4  INCIDENT  AND  THROUGHFALL 

Incident  precipitation  and  throughfall  data  for  Harp  Lake 
No.  4  have  been  compiled  in  comparable  files.  These  are: 

HL_TF_1.WK1       Intensive  Stratified  -  Reference  Nos.  36  -  103 
HL_TF_2.WK1      Reduced  Stratified  -  Reference  Nos.  104  -  186 
HL_TF_3.WK1      Plots  -  Reference  Nos.  188  -  256 

9.3  PLASTIC  LAKE  NO.  1  AND  HARP  LAKE  NO.  4  STEMFLOW 

Stemflow  was  only  collected  during  the  intensive  and  reduced 
stratified  sampling  strategies.  During  some  precipitation  peri- 
ods, there  was  no  stemflow  from  any  of  the  monitored  trees.  The 
data  files  exclude  any  reference  to  such  cases.  However,  if 
stemflow  occurred  from  at  least  one  of  the  trees,  all  trees 
being  monitored  are  listed,  frequently  with  volume  equal  to 
zero. 

Two  files  are  used  to  contain  stemflow  data.  PL_SF.WK1  con- 
tains Plastic  Lake  No.  1  data.  HL_SF.WK1  contains  Harp  Lake  No. 
4  data. 
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